The  bi -azimuthal  scattering  distribution  of  an  abyssal  hill 
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High-resolution  bistatic  images  of  a  typical  abyssal  hill  on  the  western  flank  of  the  Mid- Atlantic 
Ridge  are  made  with  a  low-frequency  towed-array  system  operating  remotely  at  1/2  convergence 
zone  (—33.3  km)  stand-off.  Comparison  with  modeled  images,  generated  from  high-resolution 
supporting  bathymetry  sampled  at  5-m  intervals,  roughly  the  wavelength  scale,  reveals  that  steep 
scarps  return  the  strongest  echoes  because  they  project  the  largest  area  along  the  acoustic  path  from 
the  source  to  receiver.  Prominent  returns  deterministically  image  scarp  morphology  when  the 
cross-range  axis  of  the  system’s  resolution  footprint  runs  along  the  scarp  axis.  Statistical  fluctuations 
inherent  in  the  scattered  field  prevent  the  system  from  distinguishing  smaller-scale  anomalies  on  the 
scarps,  such  as  canyons  and  gullies  (—100-200  m  scale),  that  would  otherwise  be  resolvable  in 
range,  in  certain  bistatic  geometries.  The  mean  bi-azimuthal  scattering  distributions  of  the  two  major 
scarps  on  the  abyssal  hill  are  identical  and  have  strengths  equal  to  the  constant  — 17  dB±8  dB.  This 
suggests  that  long-range  reverberation  from  prominent  geomorphological  features  of  the  world’s 
mid-ocean  ridges  can  be  adequately  modeled  as  Lambertian  with  albedo  7t/10L7,  given  supporting 
bathymetry  sampled  with  sufficient  frequency  to  resolve  the  projected  area  of  these  features. 
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PACS  numbers:  43.30.Hw,  43.30.Pc,  43.30.Gv  [DLB] 


INTRODUCTION 

During  the  Office  of  Naval  Research  (ONR)  Main 
Acoustics  Experiment  (MAE)  of  July  1993, 1,2  two  research 
vessels  equipped  both  with  vertical  source  and  horizontal 
receiving  arrays  made  bistatic  measurements  of  scattering 
from  a  typical,  highly  lineated,  abyssal  hill,  referred  to  as  B ' , 
on  the  western  flank  of  the  Mid- Atlantic  Ridge  (MAR).2 
These  measurements  were  made  for  two  primary  reasons. 
The  first  was  to  investigate  the  extent  to  which  an  abyssal 
hill,  one  of  the  most  commonly  occurring  bathymetric  fea¬ 
tures  of  the  world’s  oceans,  can  be  deterministically  resolved 
using  a  remotely  operated  bistatic  towed-array  imaging  sys¬ 
tem.  The  second  was  to  measure  the  bi-azimuthal  scattering 
distribution  of  a  typical  abyssal  hill. 

In  a  previous  paper,  Ref.  2,  it  was  shown  that  low- 
frequency  towed-array  systems  can  be  used  to  generate  wide- 
area  images  of  prominent  geological  features  in  the  deep 
ocean  over  hundreds  of  kilometers  in  near  real-time.  The 
brightest  features  in  these  images  effectively  resolve  the 
shape  of  steep  escarpments  (—1  km  wide)  running  along  the 
axis  of  abyssal  hills  (—10X40 km)  and  other  bathymetric 
highs  such  as  inside  corner  domes  (—10  km  diameter).2  Be¬ 
cause  the  resolution  of  these  images  greatly  exceeds  the 
roughly  10-20-km  resolution  of  currently  available  bathy¬ 
metric  maps  of  the  world’s  oceans,  where  individual  abyssal 
hills  remain  indistinguishable,  there  is  real  promise  that 
towed-array  systems  may  provide  a  rapid  means  of  charting 
previously  undiscovered  bathymetric  highs  in  the  deep 
ocean. 
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The  first  main  objective  of  this  work  is  to  experimentally 
determine  the  effective  resolving  power  of  a  towed-array 
system  in  remotely  imaging  deep-ocean  bathymetry.  To  do 
so,  an  attempt  is  made  to  infer  the  detailed  structure  of  the 
major  scarps  on  the  B '  abyssal  hill  using  MAE’s  bistatic 
system  at  1/2  convergence  zone  (CZ)  stand-off.  Since  the 
cross-range  resolution  of  the  system  (—1.0  km  @  1/2  CZ)  is 
insufficient  for  this  task,  its  much  higher  range  resolution 
(—40  m)  is  exploited.2  The  basic  idea  behind  the  analysis  is 
that  when  the  bistatic  measurement  geometry  is  such  that  the 
receiving  array’s  cross-range  direction  is  skew  to  the  ridge 
axis,  characteristic  nonlinearities  in  scarp  structure,  such  as 
canyons  (—100-200  m  wide)  and  gullies  (—50  m  wide)  pre¬ 
viously  documented  with  high-resolution  supporting 
bathymetry,2,3  become  potentially  resolvable  by  the  towed- 
array  system  in  range.  The  analysis  of  a  full  suite  of  such 
bistatic  measurements,  spanning  ±90°  about  ridge-axis 
broadside,  could  lead  to  a  detailed  accounting  of  the  scarp’s 
nonlinear  structure.  Experimental  uncertainties  in  the  orien¬ 
tation  of  the  towed-array,  sound-speed  structure  and  source- 
receiver  navigation  as  well  as  statistical  fluctuations  in  the 
scattered  returns  are  difficult  to  account  for  a  priori  and  can 
severely  degrade  the  system’s  theoretical  resolution.  The 
value  of  the  present  analysis  lies  in  its  experimental  assess¬ 
ment  of  a  towed-array  imaging  system’s  resolving  power. 

The  second  main  objective  is  to  experimentally  deter¬ 
mine  the  inherent  bi-azimuthal  scattering  distribution  of  the 
primary  scatterers  on  the  B '  abyssal  hill,  namely  two  pri¬ 
mary  scarps  that  run  along  its  major  axis,  when  waterborne 
propagation  paths  are  effectively  horizontal  at  the  abyssal 
hill,  as  they  are  in  remote  sensing  applications.  The  hope  is 
that  the  bi-azimuthal  scattering  distributions  measured  for 
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FIG.  1.  The  bistatic  tow-ship  tracks  of  the  RY  Cory  Chouest  (white  lines)  and  RV  Alliance  (black  lines)  during  the  B'  experiment  overlain  on  hydrosweep 
bathymetry  (200-m  sampling).  The  box  indicates  the  region  where  high-resolution  bathymetry  (5-m  sampling)  is  available.  Conjugate  depth  with  respect  to 
the  Cory  source  falls  at  3800  m.  Tracks  for  the  entire  B' -C'  experiment  are  shown  in  Ref.  2. 


these  may  apply  more  generally  to  scarps  found  on  abyssal 
hills  throughout  the  MAR,  many  of  which  share  the  same 
highly  lineated  character.  This  type  of  knowledge  about  the 
scattering  properties  of  abyssal  hills  is  important  because  it 
enables  one  to  (1)  efficiently  model  bistatic  scattering  from 
abyssal  hills,  which  often  appears  as  clutter  in  active  sonar 
systems  operating  in  the  deep  ocean;4  (2)  ground- truth  the 
output  of  seismo-acoustic  scattering  models  with  an  empiri¬ 
cal  scattering  function;  and  (3)  classify  bathymetric  features 
according  to  their  scattering  properties  so  similar  features 
may  be  identified  in  future  towed-array  surveys  where  there 
is  no  supporting  bathymetric  database. 

To  generalize  the  present  bistatic  scattering  measure¬ 
ments,  such  effects  as  two-way  transmission  loss,  the  areal 
resolution  of  the  measurement,  and  its  foreshortening  must 
be  factored  from  them.  In  particular,  high-resolution  bathym¬ 
etry,  sampled  at  the  scale  of  the  acoustic  wavelength,  is  used 
to  determine  surface  orientation  and  foreshortening.  The 
parabolic  equation  is  used  to  compute  two-way  transmission, 
ray  tracing  to  determine  incident  and  scattered  angles  with 
respect  to  the  seafloor  surface,  and  a  spatial  convolution  to 
account  for  the  spatially  varying  areal  resolution  of  the  mea- 
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surement  system.  Statistical  fluctuations  due  to  nonlinear 
speckle  noise  inherent  in  the  measurements  are  reduced  by 
stationary  averaging  and  are  stabilized  by  logarithmic  trans¬ 
formation.  Since  all  measurements  of  B  ’  are  at  1/2  CZ  stand¬ 
off,  acoustic  paths  from  the  source  and  to  the  receiver  are 
nearly  horizontal  when  they  intersect  B ' .  This  enables  a  rela¬ 
tively  accurate  calculation  of  the  bi-azimuthal  scattering  dis¬ 
tribution  of  the  B  ’  scarps  which  will  be  useful  in  the  classi¬ 
fication  of  similar  scarps  on  mid-ocean  abyssal  hills  by 
towed-array  systems  at  n  +  1/2  CZ. 

I.  BISTATIC  EXPERIMENTAL  DESIGN,  IMAGING 
RESOLUTION,  AND  THE  GEOMORPHOLOGY  OF  THE 
B'  ABYSSAL  HILL 

The  geometry  of  the  bistatic  experimental  survey  of  the 
B'  abyssal  hill  is  summarized  here  in  Fig.  1,  where  the  track 
design,5,6  is  overlain  on  hydrosweep  bathymetry,  sampled  at 
200-m  intervals.2,7  A  more  complete  description  of  the 
B’  -C'  corridor  experiments,  which  comprised  roughly  90% 
of  MAE,  can  be  found  in  Refs.  2,  5,  and  6.  The  8X9-km2 
box  at  the  center  of  B'  shows  the  region  where  high- 
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FIG.  2.  (a)  High-resolution  bathymetry  (5-m  sampling)  on  the  east-central  face  of  B'  as  located  in  Fig.  1.  Depth  contours  are  at  100-m  intervals,  (b)  The 
component  of  seafloor  elevation  gradient,  or  the  directional  derivative  (DD),  normal  to  the  B'  axis,  computed  from  high-resolution  bathymetry  as  in  Fig.  22 
of  Ref.  2.  Many  steep  cliffs  (>60°  slope)  and  smaller  terraces  (—45°  slope)  appear  on  both  the  upper  and  lower  scarps  that  cannot  be  resolved  by  the 
hydrosweep  bathymetry.  Contours  segment  the  upper  and  lower  scarps  as  regions  within  which  the  DD  exceeds  1/2.  These  same  contours  are  used  to  identify 
the  scarps  throughout  this  paper. 


resolution  bathymetry,  sampled  at  5-m  intervals,  was  col¬ 
lected  in  support  of  MAE.2,3  The  research  vessel  (RV)  CORY 
Chouest  traced  the  central  star  while  the  RV  Alliance 
followed  the  wings  to  the  north  and  south.  Each  RV  towed  a 
horizontal  receiving  array  with  axis  roughly  coincident  with 
the  ship’s  straight-line  course.  To  maximize  cross-range 
resolution  of  B'  and  to  minimize  noise  contamination  from 
the  other  RV,  the  highest  resolution  beams  at  array  broadside 
are  directed  towards  the  center  of  B ' ,  where  high-resolution 
bathymetry  is  available,  while  the  lowest  resolution  beams  at 
array  endfire  are  directed  towards  the  other  RV. 2,4,6  Only 
transmissions  from  the  CORY  CHOUEST  source  array  are  con¬ 
sidered  in  this  paper  because  of  its  significantly  higher 
strength,  229  dB  re  1  mPa  @  1  m,  and  directionality  com¬ 
pared  to  that  of  the  Alliance.2,8  The  Cory’s  source  array  is 
steered  to  broadside  for  all  transmissions  studied  in  this  pa¬ 
per. 

Linear  frequency  modulated  (LFM)  waveforms  of  r 
=  5  s  duration  in  the  frequency  band  200-255  Hz  are  exclu¬ 
sively  analyzed  here  because  they  offer  the  highest  available 
range  resolution,  which  is  given  by  Ar  =  c/(2B)~  14m, 
where  15 10  m/s  is  the  mean  sound  speed  and  B  is  the 
bandwidth.2  To  reduce  their  standard  deviation  after  pulse 
compression,  instantaneous  intensity  measurements  are  aver¬ 
aged  over  time  period  T  =  0.0625  s  for  CORY  receptions  and 
T=  0.0533  s  for  Alliance  receptions.  This  averaging  leads 
to  an  effective  range  resolution  of  Ar  =  c  772 =47  m  for  the 
Cory  and  A  r  =  40  m  for  the  Alliance. 

With  the  assumption  that  the  specific  bathymetric  fea¬ 
tures  resolved  by  the  towed-array  system  behave  as  Rayleigh 
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targets,9,10  due  to  their  structural  complexity  with  respect  to 
the  wavelength  scale,  the  degrees-of-freedom  N~TB  be¬ 
comes  roughly  3.4  for  Cory  and  2.9  for  Alliance  LMF 
receptions.  The  further  assumption  of  independence  and  sta- 
tionarity  among  the  roughly  3  range  cells  averaged  leads  to  a 
standard  deviation  of  2.5  dB  and  a  bias  of  0.6  dB  in  the 
intensity  level  at  any  pixel  of  a  towed-array  image.2,10  If,  on 
the  other  hand,  a  particular  range  cell  dominates  the  average, 
or  if  the  range  cells  are  correlated,  /z  may  be  as  small  as 
unity,  which  leads  to  a  standard  deviation  of  roughly  5.6  dB 
and  a  bias  of  2.5  dB.2,10 

The  geomorphology  of  the  B'  abyssal  hill  is  displayed 
in  Fig.  2.  Gross  features  are  evident  in  the  high-resolution 
bathymetry  of  Fig.  2(a).  The  central  crest  of  B ' ,  for  example, 
rises  roughly  1400  m  from  the  western  segment  valley  floor 
at  roughly  4900  m  to  nearly  400  m  above  the  3800-m  con¬ 
jugate  depth2  of  the  CORY  source.  Minute  details  of  the  high- 
resolution  bathymetry  are  best  cast  in  terms  of  the  compo¬ 
nent  of  seafloor  elevation  gradient,  or  directional  derivative, 
normal  to  the  B'  axis,  shown  in  Fig.  2(b).  This  reveals  two 
steep  scarps  separated  by  an  extended  plateau  just  below 
3800  m.  Both  scarps  run  along  the  B'  axis  21°  from  true 
North,  and  are  characterized  by  steep  unsedimented  cliff 
faces  that  rise  upwards  for  typically  100  m  with  characteris¬ 
tic  slopes  of  50° -90°.  The  scarps  are  segmented  by  a  contour 
within  which  the  directional  derivative  exceeds  1/2,  corre¬ 
sponding  to  slope  angles  greater  than  —26.6°.  The  contours, 
shown  in  Fig.  2(b),  are  used  throughout  the  paper  to  desig¬ 
nate  the  locations  of  the  upper  and  lower  scarps. 

The  bistatic  locations  of  the  two  research  vessels  during 
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FIG.  3.  The  bistatic  locations  of  the  two  research  vessels  during  the  LFM 
transmissions  analyzed  in  this  paper  are  given  in  eastings,  x,  and  northings, 
y.  These  locations  are  distributed  in  a  semicircle  about  the  high-resolution 
region  shown  by  the  central  box  within  which  two  diagonal  lines  indicate 
the  upper  and  lower  scarps.  Circled  letters  indicate  the  location  of  the  Cory 
while  the  same  letters  indicate  the  location  of  the  Alliance  for  the  20 
bistatic  segments  given  in  the  legend.  Segments  S430,  S466,  S481,  and 
S495  are  excluded  from  the  monostatic  analysis  due  to  colocation  with  some 
of  the  other  monostatic  segments  shown.  Since  Alliance  was  at  the  apex  of 
the  northern  wing  and  in  the  shadow  of  B'  during  S445  and  S450,  they  are 
used  as  the  only  purely  monostatic  transmissions  in  this  analysis,  which  then 
has  a  total  of  18  monostatic  segments. 


the  LFM  transmissions  analyzed  in  this  paper  are  given  in 
eastings  and  northings  in  Fig.  3.  These  locations  are  distrib¬ 
uted  in  a  semicircle  with  bisecting  diameter  oriented  along 
the  B'  axis  and  centered  in  the  high-resolution  region.  A 
letter  of  the  alphabet  designates  each  transmission.  Circled 
letters  indicate  the  monostatic  location  of  the  CORY  while 
unboxed  letters  indicate  the  bistatic  location  of  the  Alliance 
for  a  given  transmission.  A  segment  number  specifies  each 
transmission,  following  the  vernacular  of  the  MAE.1,2 

The  portion  of  bi-azimuthal  space  spanned  by  the 
source -receiver  pairs  is  shown  in  Fig.  4.  Azimuth  is  mea¬ 
sured  counterclockwise  from  a  vector  that  originates  from 
the  center  of  the  given  scarp  in  the  high-resolution  region 
and  points  normal  to  its  axis  in  a  southeasterly  direction.  The 
source  azimuths  12 ;  then  fall  within  |I2;|<30°,  while  the 
receiver  azimuths  II r  span  the  full  180°  range.  This  bistatic 
geometry  is  appropriate  because  (1)  the  east-central  face  of 
B'  is  uniformly  insonified  with  negligible  shadowing,  and 
(2)  the  receivers  are  never  completely  shadowed. 

Expressions  for  the  half-power  or  3-dB  beamwidth  (5(6) 
of  a  narrow-band  line  array,  steered  to  an  angle  6  from  end- 
fire,  are  given  in  Eqs.  (1)  and  (2)  of  Ref.  11  as  a  function  of 
aperture  length  and  frequency.  These  expressions  also  ap¬ 
proximate  the  azimuthal  resolution  of  the  MAE  arrays  for 
LFM  reception  when  evaluated  at  the  LFM  center  frequency 
of  227.5  Hz.  Best  resolution  lies  at  broadside  and  is  —1.6° 
for  the  L  =  318  m  receiving  array  of  the  CORY  and  —2.0°  for 
the  L  =  254  m  receiving  array  of  the  Alliance.  Cross-range 
resolution,  defined  as  r(5(0)  for  a  monostatic  measurement 
at  range  r  from  the  scattering  patch,  is  then  0.9  km  for  the 


FIG.  4.  The  portion  of  bi-azimuthal  space  spanned  by  the  source-receiver 
pairs.  Azimuth  is  measured  counterclockwise  from  a  vector  that  originates 
from  the  center  of  the  given  scarp  in  the  high-resolution  region  and  points 
normal  to  its  axis  in  a  southeasterly  direction.  Source  azimuths  are  denoted 
by  11 1  and  receiver  azimuths  by  Hr .  The  curves  for  the  upper  and  lower 
scarps  are  nearly  identical  because  their  separation  is  small  compared  to  the 
range  to  the  respective  RV.  The  portion  of  nonmonostatic  space  spanned  is 
double  that  shown  by  acoustic  reciprocity  between  source  and  receiver  lo¬ 
cations. 


Cory  and  1.1  km  for  the  Alliance  at  1/2  CZ  where  r 
=  33.3  km. 

At  1/2  CZ,  neither  array  has  sufficient  aperture  to  azi- 
muthally  resolve  characteristic  nonlinearities  in  scarp  struc¬ 
ture,  such  as  canyons  (—100-200  m  wide)  and  gullies  (—50 
m  wide)  previously  documented  with  high-resolution  sup¬ 
porting  bathymetry.2,3  This  is  evident  upon  inspection  of  Fig. 
5(a),  where  the  resolution  footprint  of  the  towed- array  sys¬ 
tem,  a  roughly  50  m  by  1  km  annular  sector  at  1/2  CZ,  is 
superimposed  on  a  bathymetric  image  of  the  B '  upper  scarp. 
These  same  features  become  potentially  resolvable  by  the 
towed-array  system  in  range  when  cross-range  is  skew  to  the 
ridge  axis,  as  shown  in  Fig.  5(b),  with  the  relatively  safe 
assumption  that  the  scarp,  and  not  the  plateaus  above  and 
below,  cause  the  predominant  reverberation.  Analysis  of  a 
full  suite  of  such  bistatic  measurements,  spanning  |ftr| 
<90°,  could  provide  a  detailed  accounting  of  the  scarp’s 
nonlinear  structure,  given  sufficient  charting  accuracy  and 
statistical  stability  of  the  measured  reverberation. 

The  normalized  autocorrelation  function  of  upper  scarp 
elevation,  in  Fig.  6,  quantitatively  reveals  along  and  across- 
scarp  correlation  scales  (e-folding  lengths)  of  150  and  50-m, 
respectively,  consistent  with  a  visual  assessment  of  typical 
canyon  and  gully  dimensions.  The  resolution  footprint  of  the 
towed-array  system  integrates  or  blurs  together  scattered  re¬ 
turns  from  roughly  five  major  canyons  when  its  cross-range 
range  axis  parallels  the  scarp’s,  as  is  characteristic  of  the 
monostatic  measurements.  This  leads  to  a  broad  distribution 
of  incident  angles  from  a  given  source,  or  scattered  angles  to 
a  given  receiver,  over  the  resolution  footprint.  This  is  shown 
in  Fig.  7(a)  for  S435  where  vertical  refraction  is  taken  into 
account  via  raytrace.  Conversely,  the  resolution  footprint 
typically  blurs  together  scattered  returns  from  the  outer  wall 
of  a  single  canyon  with  those  from  plateaus  above  and  below 
when  its  cross-range  axis  is  skew  to  the  scarp  axis,  as  is 
characteristic  of  bistatic  measurements  with  large  |I2r|.  This 
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FIG.  5.  Shaded  relief  of  a  portion  of  the  B'  upper  scarp  imaged  with  high- 
resolution  bathymetry  with  the  roughly  50-m  by  1-km  resolution  footprint  of 
the  towed- array  system  overlain,  (a)  Monostatic  geometry  for  S435.  (b) 
Bistatic  geometry  for  S435.  In  (a),  the  resolution  footprint  blurs  together 
many  characteristic  nonlinearities  in  scarp  structure,  such  as  canyons  and 
gullies.  These  features  become  potentially  resolvable  in  range,  as  in  (b),  for 
some  bistatic  geometries  of  the  present  experiment. 


also  leads  to  a  broad  distribution  of  incident  angles  or  scat¬ 
tered  angles  over  the  resolution  footprint,  as  shown  in  Fig. 
7(b)  and  (c)  for  the  same  segment. 

Since  the  surfaces  resolved  by  the  towed-array  system 
on  the  B  ’  scarps  are  not  even  approximately  planar,  there  is 
no  unique  surface  normal  to  characterize  the  bathymetry 
within  the  resolution  footprint.  Consequently,  there  is  no 
unique  incident  or  scattered  angle.  The  traditional  surface 
scattering  strength ,  which  depends  upon  the  directions  of 
incident  and  scattered  waves  defined  relative  to  the  local  sur¬ 
face  normal  and  includes  purely  local  foreshortening  effects , 
must  then  be  reanalyzed  to  accommodate  the  reality  of  the 
present  experiment,  as  will  be  discussed  in  later  sections. 
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FIG.  6.  The  normalized  2-D  autocorrelation  of  upper  scarp  elevation  as  a 
function  of  east  and  north  spatial  lag.  The  along-scarp  and  across-scarp 
e-folding  lengths  of  —150  and  —50  m,  respectively,  quantify  the  scale  of 
typical  nonlinearities  along  and  across  the  scarp. 


S478  and  S487.  Analysis  of  S478  is  particularly  significant 
because  it  characterizes  all  measurements  for  which  the  bi¬ 
static  separation  between  the  Cory  and  Alliance  is  small 
relative  to  the  range  to  B  ’ .  In  these  cases,  bistatic  and  mono¬ 
static  charts  are  similar  with  cross-range  falling  parallel  to 
the  B'  axis.  Prominent  reverberation  is  then  unambiguously 
charted  to  the  steep  escarpments,  as  demonstrated  in  Ref.  2. 
Monostatic  and  bistatic  reverberation  charts  of  the  same 
transmission  show  great  differences,  however,  when  the  bi¬ 
static  separation  is  large.  This  is  well  illustrated  here  for 
S435  and  S528,  which  are  respectively  at  the  northern  ex¬ 
treme  and  midway  to  the  southern  extreme  of  the  Alliance 
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II.  BISTATIC  IMAGING  OF  THE  B'  ABYSSAL  HILL 
FROM  1/2  CZ 

A.  Wide-area  bistatic  images 

While  a  number  of  monostatic ,  reverberation  charts  of 
the  B'  abyssal  hill  have  been  presented  and  analyzed, 2,1 2-1 4 
only  two  bistatic  charts  of  B '  have  been  previously  studied 
in  the  literature.2  These  are  for  the  200-255-Hz  LFMs  of 


FIG.  7.  The  distribution  of  incident  and  scattered  angles  within  the  resolu¬ 
tion  footprints  shown  in  Fig.  5  for  S435.  The  incident  angle  is  the  inverse 
cosine  of  the  inner  product  between  the  local  surface  normal  and  the  vector 
pointing  along  the  acoustic  path  to  the  source.  The  scattered  angle  is  the 
inverse  cosine  of  the  inner  product  between  the  local  surface  normal  and  the 
vector  pointing  along  the  acoustic  path  to  the  receiver.  Refraction  is  in¬ 
cluded.  Angles  greater  than  90°  indicate  shadowing,  (a)  Monostatic  case 
where  incident  and  scattered  angles  are  identical  and  axis  of  resolution 
footprint  is  along  scarp  axis.  Bistatic  case  for  (b)  incident  angles  and  (c) 
scattered  angles,  where  axis  of  resolution  footprint  is  skew  to  scarp  axis. 
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FIG.  8.  Wide-area  charts  of  mono¬ 
static  and  bistatic  reverberation  mea¬ 
sured  for  the  200-255-Hz  LFM  S435, 
where  Alliance  is  near  the  apex  of 
the  northern  wing,  (a)  Monostatic  re¬ 
verberation  chart  showing  symmetry 
about  the  array  axis  for  Cory  heading 
163°.  (b)  Contours  of  high-level  back- 
scatter,  overlain  on  the  bistatic  hori¬ 
zontal  projection  of  bathymetry 
(BHPB),  coregister  with  major  scarps 
on  B'  that  have  large  areas  projected 
towards  the  source-receiver.  The  false 
backscatter  image  of  B '  does  not 
coregister  with  the  BHPB.  (c)  Bistatic 
reverberation  chart  shows  asymmetry 
about  the  array  axis  because  this  axis 
forms  an  oblique  angle  with  the  major 
axis  of  the  bistatic  ellipse.  The  Alli¬ 
ance  heading  is  277°.  (d)  Contours  of 
high-level  backscatter,  overlain  on  the 
BHPB,  coregister  with  major  scarps 
on  B'  that  have  large  areas  projected 
along  the  acoustic  path  from  source  to 
scatterer  to  receiver.  Unlike  the  mono¬ 
static  case,  the  cross-range  axis  of  the 
resolution  footprint  forms  an  oblique 
angle  with  the  scarp  axis  so  that 
prominent  returns  are  speckled  across 
the  scarp.  The  false  backscatter  image 
of  B'  does  not  coregister  with  the 
BHPB  and  exhibits  severe  spatial  con¬ 
traction.  Source  and  receiver  locations 
are  indicated  in  (a)  and  (c).  The  prod¬ 
uct  of  a  negative  HPBR  and  a  negative 
HBPS  is  always  plotted  as  a  negative 
HPBR  to  distinguish  shadowed  re¬ 
gions. 


tracks.  Specifically,  - 5°  and  flr~90°  for  S435  and 

fl, - 5°  and  £lr - 50°  for  S528,  with  respect  to  both  up¬ 

per  and  lower  scarp  reference  centers. 

To  gain  a  broader  perspective,  consider  first  the  wide- 
area  charts  of  monostatic  reverberation,  in  Figs.  8- 10(a). 
These  each  have  a  character  very  similar  to  that  found  for 
short  continuous  wave  (cw)  transmissions  in  Figs.  8  and  9  of 
Ref.  2.  When  overlain  on  the  bistatic  horizontal  projection  of 
bathymetry  (BHPB),  prominent  reverberation  coregisters 
precisely  with  steep  escarpments  along  the  B '  axis,  as  illus¬ 
trated  in  Figs.  8- 10(b).  The  bistatic  horizontal  projection  of 
bathymetry  is  the  product  of  the  horizontal  projection  of 
bathymetry  toward  the  source  (HPBS)  and  the  horizontal 
projection  of  bathymetry  toward  the  receiver  (HPBR).  The 
former  is  inner  product  of  the  local  seafloor  normal  and  a 
unit  vector  pointing  to  the  horizontal  position  of  the  source, 
while  the  latter  is  the  inner  product  of  the  local  seafloor 
normal  and  a  unit  vector  pointing  to  the  horizontal  position 
of  the  receiver.  To  distinguish  shadowed  regions,  the  product 
of  a  negative  HPBR  and  a  negative  HBPS  is  always  plotted 
as  a  negative  HPBR.  In  later  sections,  a  far  more  precise 
computation  of  surface  projection  factors,  accounting  for  re¬ 
fracted  ray  paths  to  and  from  the  seafloor,  is  employed  to 
analyze  reverberation  over  the  region  of  high-resolution 
bathymetry. 

In  the  general  limiting  case  of  range-independent  propa- 
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gation  and  a  flat  bottom,  range  ambiguity  of  the  bistatic 
towed-array  system  falls  on  elliptic  arcs  about  foci  coincid¬ 
ing  with  the  source  and  receiver  positions.  The  arcs  are  de¬ 
fined  by  conserving  the  two-way  travel  distance  from  source 
to  scatterer  to  receiver  as  the  sum  of  the  radii.  The  arc  length 
at  1/2  CZ,  over  the  broadside  resolution  of  the  CORY  or 
Alliance  receiving  array,  remains  well  approximated  by  the 
monostatic  cross-range  resolution  rfi(0)  for  the  entire  suite 
of  bistatic  measurements  taken  at  B ' ,  as  shown  in  Appendix 
A.  Here,  r  is  range  from  the  receiving  array  to  the  scattering 
site.  For  what  we  refer  to  here  as  ‘‘monostatic”  measure¬ 
ments,  the  Cory’s  source  and  receiver  arrays  are  not  coinci¬ 
dent  but  are  in  such  close  proximity  (~1.12-km  separation) 
that  range  ambiguity  lies  on  nearly  circular  arcs  at  1/2  CZ. 
When  the  bistatic  separation  between  the  CORY  and  ALLI¬ 
ANCE  is  not  small  compared  to  r,  however,  the  arcs  are  far 
from  circular.  This  accounts  for  some  characteristic  differ¬ 
ences  between  monostatic  and  bistatic  charts  of  the  same 
transmission. 

The  right-left  ambiguity  of  the  receiving  line  array  is 
also  expressed  differently  in  monostatic  and  bistatic  charts. 
Prominent  features  are  ambiguously  charted  in  an  almost 
symmetrical  fashion  about  the  receiving  array’s  axis  in 
monostatic  geometries.11,2  Only  the  incorporation  of  known 
bathymetry  in  the  charting  procedure  leads  to  local  breaks  in 
symmetry.2,14  Such  extreme  symmetry,  however,  is  only 
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FIG.  9.  Wide-area  charts  of  mono¬ 
static  and  bistatic  reverberation  mea¬ 
sured  for  the  200-255-Hz  LFM  S478. 
Alliance  and  Cory  are  in  close  prox¬ 
imity.  (a)  Monostatic  reverberation 
chart  showing  symmetry  about  the  ar¬ 
ray  axis  for  Cory  heading  45°.  (b) 
Contours  of  high-level  backscatter, 
overlain  on  the  BHPB,  coregister  with 
major  scarps  on  B '  that  have  large  ar¬ 
eas  projected  towards  the  source- 
receiver.  (c)  Bistatic  reverberation 
chart  shows  symmetry  about  the  array 
axis  due  to  the  close  proximity  of  the 
Cory  and  Alliance.  The  Alliance 
heading  is  207°.  (d)  Contours  of  high- 
level  backscatter,  overlain  on  the 
BHPB,  coregister  with  major  scarps 
on  B'  that  have  large  areas  projected 
along  the  acoustic  path  from  source  to 
scatterer  to  receiver. 


available  in  bistatic  scenarios  when  the  receiving  array  axis 
parallels  the  line  joining  the  source  and  receiver.  This  line 
defines  the  major  axis  of  the  bistatic  ellipse.  Otherwise,  some 
distortion  occurs.  In  the  case  of  segment  435,  the  receiving 
array  is  at  an  oblique  angle  to  the  major  axis  of  the  bistatic 
ellipse.  Prominent  reverberation  from  the  central  scarps  of 
B'  correctly  register  at  roughly  1/2  CZ  (—33.3  km)  from  the 
receiving  array.  The  corresponding  ambiguous  returns  are 
mirrored  about  the  array  axis,  and  span  the  same  annular 
sector  as  the  true  returns,  but  fall  at  a  much  shorter  range  to 
preserve  the  travel  time.  As  a  result,  they  are  condensed  to  a 
much  smaller  spatial  area  than  that  spanned  by  the  true  re¬ 
turns.  The  converse  is  also  true.  True  returns  from  this  small 
area,  where  the  right-left  ambiguity  of  the  receiving  array 
intersects  the  elliptical  range  ambiguity  of  the  bistatic  sys¬ 
tem,  may  be  falsely  charted  to  a  more  extensive  region  on 
B ' .  The  geometry  of  this  type  of  ambiguity  is  given  in  Fig. 
Al.  After  a  careful  analysis  of  the  registration  between 
bathymetric  directional  derivatives  and  prominent  reverbera¬ 
tion,  we  find  that  the  experimental  design5  was  successful  in 
exploiting  the  natural  geomorphology  of  the  region  to  avoid 
contamination  from  line-array  ambiguity  at  B ' . 

Prominent  bistatic  returns  are  overlain  on  the  bistatic 
horizontal  projection  of  bathymetry  (BHPB)  in  Figs. 
8- 10(d).  In  the  limit  of  horizontal  propagation  to  and  from 
the  bottom,  which  is  a  reasonable  approximation  for  water- 
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borne  paths  at  n  + 1/2  CZ,  the  BHPB  becomes  a  projection 
factor  directly  proportional  to  the  incoherent  intensity  of  the 
scattered  field  from  an  infinitesimal  area  of  the  seafloor,  and 
so  is  useful  for  comparison  with  wide-area  reverberation. 

For  large  bistatic  separations,  as  in  S435  and  S528,  char¬ 
acteristic  nonlinearities  in  scarp  structure,  such  as  canyons 
and  gullies,  lead  to  a  speckled  pattern  of  strong  projections  in 
the  path  connecting  source,  scatterer,  and  receiver.  Over  the 
same  region,  prominent  bistatic  reverberation  displays  analo¬ 
gous  behavior  as  a  speckled  pattern  of  concentric  arcs  that 
obliquely  cross  the  scarp  axes.  This  phenomenon  was  previ¬ 
ously  observed  in  monostatic  reverberation  charts  when  the 
radial  propagation  path  was  parallel  to  or  formed  a  shallow 
angle  with  the  axis  of  a  major  bottom-limited  ridge.2  The 
effect  is  partially  due  to  (1)  anomalous  surface  projections, 
and  (2)  statistical  fluctuations  from  random  interference,  as 
will  be  shown  in  the  next  sections. 

B.  High-resolution  images  of  the  east-central  face  of 

B' 

The  bending  of  acoustic  paths  from  the  source,  to  the 
scattering  patch  and  back  to  the  receiver,  caused  by  naturally 
occurring  depth-dependent  variations  in  the  ocean’s  refrac¬ 
tive  index,  described  in  Ref.  2,  must  be  taken  into  account 
for  the  charting  accuracy  of  long-range  echoes  to  be  compa- 
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FIG.  10.  Wide-area  charts  of  mono¬ 
static  and  bistatic  reverberation  mea¬ 
sured  for  the  200-255-Hz  LFM  S528. 
Alliance  is  midway  along  the  southern 
wing,  (a)  Monostatic  reverberation 
chart  showing  symmetry  about  the  ar¬ 
ray  axis  for  Cory  heading  206°.  (b) 
Contours  of  high-level  backscatter, 
overlain  on  the  BHPB,  coregister  with 
major  scarps  on  B '  that  have  large  ar¬ 
eas  projected  towards  the  source- 
receiver.  The  false  backscatter  image 
of  B'  does  not  coregister  with  the 
BHPB.  (c)  Bistatic  reverberation  chart 
shows  asymmetry  about  the  array  axis 
due  to  the  oblique  angle  it  forms  with 
the  major  axis  of  the  bistatic  ellipse. 
The  Alliance  heading  is  87°.  (d) 
Contours  of  high-level  backscatter, 
overlain  on  the  BHPB,  coregister  with 
major  scarps  on  B '  that  have  large  ar¬ 
eas  projected  along  the  acoustic  path 
from  source  to  scatterer  to  receiver. 
The  false  backscatter  image  of  B' 
does  not  coregister  with  the  BHPB  and 
exhibits  severe  spatial  contraction. 


rable  to  the  5-m  sample  interval  of  the  high-resolution 
bathymetry  shown  in  Fig.  2.  A  sophisticated,  range- 
dependent  ray-based  propagation  model  known  as  RASP,15 
that  includes  the  effects  of  multiple  interactions  between  sur¬ 
face  and  bottom  boundaries,  was  used  to  chart  reverberation 
out  to  multiple  CZ  ranges  in  Ref.  2.  A  number  of  other 
methods  have  been  developed  and  employed  by  other  SRP 
researchers.13,14,16 

For  the  purpose  of  the  present  analysis,  where  only 
waterborne  insonification  of  the  high-resolution  area  is  con¬ 
sidered,  a  somewhat  simpler  but  more  efficient  raytrace 
model  is  used  to  obtain  the  two-way  travel  times.  Rays  are 
traced  in  a  range-depth  plane  from  a  point  source,  located  at 
the  center  of  the  respective  source  or  receiving  array,  through 
a  refractive  ocean  half-space  with  a  pressure-release  surface 
but  no  bottom  boundary.  Surface-reflected  paths  are  ac¬ 
counted  for,  but  the  much  slower  and  more  attenuated 
bottom-reflected  paths  are  neglected.  Range-depth  maps  of 
minimum  travel  time  are  then  generated.  A  range-depth  map 
is  then  swept  through  the  high  resolution  bathymetry  of  B  ’ 
to  determine  minimum  travel  time  to  or  from  the  bottom. 
This  procedure  is  repeated  for  all  source  and  receiver  loca¬ 
tions  to  obtain  respective  travel  time  maps  in  the  horizontal 
plane  over  the  high-resolution  region.  These  maps  are  then 
used  to  chart  the  output  of  the  beamformed  and  matched 
filtered  reverberation  data  to  the  high-resolution  region  by 
the  same  bistatic  mapping  algorithm  employed  in  Refs.  2, 
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11,  and  12.  Navigation  data  from  the  Global  Positioning  Sys¬ 
tem  (GPS)17  of  each  RV  are  precisely  converted  to  universal 
transverse  Mercator  (UTM)18  coordinates,  for  long-range  re¬ 
verberation  charting,  by  the  United  States  Geological  Survey 
(USGS)  geodetic  reference  system19  that  is  accurate  to 
within  roughly  1  m. 

The  error  in  computing  travel  time  to  the  high-resolution 
region  by  slant-range,  rather  than  raytrace,  is  illustrated  in 
Fig.  11.  Since  the  one-way  errors  shown  are  typically  much 
greater  than  the  roughly  0.003  s  it  takes  to  cross  a  5-m 
sample  interval  of  high-resolution  bathymetry  at  the  mean 
sound  speed,  the  slant-range  approach  to  charting  reverbera¬ 
tion  is  unacceptable. 

High-resolution  images  of  charted  reverberation  are  il¬ 
lustrated  in  Figs.  12-14  for  segments  435,  478  and  528,  with 
contours  overlain  to  specify  locations  of  the  upper  and  lower 
scarps.  Both  scarps  tend  to  appear  prominently  and  faithfully 
imaged  in  the  monostatic  charts,  as  in  S478  and  S435.  This 
fidelity  is  possible  because  the  cross-range  axis  of  the  reso¬ 
lution  footprint  lies  along  the  scarp  axes.  In  some  segments, 
such  as  S528,  however,  the  entire  upper  scarp  vanishes 
amidst  low-level  background  reverberation.  This  happens 
when  the  upper  scarp  falls  in  the  shadow  zone  of  the  main 
beam’s  refractive  path,  due  to  its  proximity  to  the  source 
array’s  conjugate  depth.2  The  lower  scarp  appears  promi¬ 
nently  in  all  monostatic  charts  because  it  lies  below  the 
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FIG.  11.  (a)  Slant-range  travel  time,  using  mean  sound  speed  along  the  slant  path,  minus  raytrace  travel  time  from  Cory  source  to  seafloor  for  S435  using 
high-resolution  bathymetry,  (b)  Same  as  (a)  except  from  seafloor  to  Alliance  receiver  for  S435.  These  one-way  errors  are  typically  much  greater  than  the 
0.003-s  travel  time  it  takes  to  cross  a  5-m  sample  interval  of  high-resolution  bathymetry. 


shadow  zone,  displayed  in  Fig.  15(a),  and  generally  receives 
main-beam  insonification. 

The  consistent  coregistration  between  the  locations  of 
the  strongest  monostatic  returns  and  the  locations  of  the 
scarps  carries  much  statistical  weight.  The  upper  5%  of  re¬ 
verberation  within  the  high-resolution  area  is  selected  for 
each  segment  to  demonstrate  this.  The  frequency  of  overlap 
of  these  selected  regions,  over  all  18  monostatic  segments,  is 
plotted  in  Fig.  16(a)  together  with  overlays  of  the  upper  and 
lower  scarps.  This  figure  shows  that  the  regions  that  most 
frequently  return  the  strongest  backscatter  coincide  with 
some  of  the  steepest  portions  of  the  upper  and  lower  scarps. 
The  figure  also  shows  that  the  frequency  of  strong  backscat¬ 
ter  rapidly  decreases  as  the  distance  from  the  scarp  increases 
in  the  cross-scarp  direction.  It  is  extremely  unlikely,  for  ex¬ 
ample,  to  find  even  a  single  strong  return  charted  to  points 
lying  at  ranges  beyond  roughly  1.0  km  of  either  scarp. 

The  bistatic  charts  for  S435  and  S528  are  significantly 
different  from  the  corresponding  monostatic  charts,  as  noted 
in  the  wide-area  analysis  of  the  previous  section.  The  differ¬ 
ence  arises  because  the  system’s  cross-range  resolution  is  at 
an  oblique  angle  to  the  scarp  axes  in  the  bistatic  case, 
whereas  it  parallels  it  in  the  monostatic  case.  Prominent  bi¬ 
static  returns  lie  along  nearly  elliptical  arcs  centered  across 
the  scarp  axes.  This  is  seen  in  Figs.  12-14,  which  are  typical 
of  many  of  the  bistatic  segments  studied.  The  strongest  re¬ 
turns  appear  in  a  speckled  fashion  along  the  scarps  axis.  In 
extreme  cases  where  cross-range  is  nearly  normal  to  the 
scarps  axes,  as  in  S435,  the  towed- array  system  is  sometimes 
unable  to  distinguish  the  upper  from  the  lower  scarp.  This 
effect  is  apparently  due  to  the  relatively  close  proximity  of 


the  two  scarps,  roughly  2  km  separation  or  twice  the  cross¬ 
range  resolution  of  the  receiver,  and  occasionally  leads  to 
strong  returns  continuously  charted  across  both  scarps  in  el¬ 
liptic  arcs  centered  through  the  plateau.  This  may  lead  one  to 
the  presumably  false  conclusion  that  the  plateau  contains 
geomorphological  anomalies  that  may  yield  strong  returns. 

The  spatial  frequency  distribution  of  the  most  prominent 
bistatic  returns  is  shown  in  Fig.  16(c).  As  in  the  monostatic 
case,  the  upper  5%  of  reverberation  over  the  high-resolution 
area  is  selected  for  each  segment.  The  overlap  frequency  of 
these  selected  regions,  across  all  20  bistatic  segments,  is 
plotted  together  with  overlays  of  the  upper  and  lower  scarps. 
While  the  global  maxima  appears  on  the  steepest  portion  of 
the  upper  scarp,  and  local  maxima  often  register  with  both 
the  upper  and  lower  scarps,  the  largest  number  of  local 
maxima  appear  in  the  region  between  the  scarps.  Frequency 
tends  to  decrease  as  range  increases  from  the  region  between 
and  including  both  scarps.  Single-frequency  returns,  how¬ 
ever,  remain  densely  spread  with  relative  uniformity  from 
the  upper  crest  to  the  lower  scarp,  and  only  begin  to  disap¬ 
pear  as  the  trench  below  the  lower  scarp  is  approached. 
While  some  cross-hatching  of  overlain  returns  is  to  be  ex¬ 
pected  just  above  and  below  the  upper  and  lower  scarps,  due 
to  a  combination  of  ridge  morphology,  track  geometry,  and 
cross-range  blurring,  the  significantly  higher  concentration  of 
high  frequencies  between  the  upper  and  lower  scarps  rules 
this  combined  effect  out  as  the  sole  cause  of  off- scarp  chart¬ 
ing  of  prominent  reverberation,  as  is  shown  in  the  next  sec¬ 
tion.  An  explanation  for  the  measured  frequency  distribution 
follows  in  Sec.  II D. 
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FIG.  12.  High-resolution  charts  of 
monostatic  and  bistatic  measured  and 
modeled  reverberation  for  S435  over 
the  region  shown  in  Fig.  2  with  upper 
and  lower  scarp  contours  overlain,  (a) 
Measured  monostatic  reverberation. 

(b)  Modeled  monostatic  reverberation. 

(c)  Measured  bistatic  reverberation. 

(d)  Modeled  bistatic  reverberation. 


C.  Modeling  high-resolution  reverberation  and 
estimating  the  scattering  distribution 

From  an  infinitesimal  planar  surface  patch  of  area  dA ,  in 
the  far  field  of  a  source  and  receiver,  the  scattered  intensity  at 
the  receiver  can  be  written  as 

dl= f(0i  ,4>r)wljl r  cos  6 1  cos  6rdA ,  (1) 

according  to  the  basic  principles  of  radiometry.  Here,  0  and 
4>  denote  polar  and  azimuth  angles  with  respect  to  the  local 
surface  normal,20  with  incident  and  scattered  angles  identi¬ 
fied  by  the  subscripts  i  and  r.  The  projected  area  of  the  sur¬ 
face  patch  in  the  direction  of  the  scattered  wave  leaving  the 
surface  patch  is  cos  6rdA.  The  other  factors  include  the 
source  power  w,  transmission  coefficients  lt  and  lr  in 
inverse- square  meters  from  the  source  to  the  surface  patch 
and  from  the  surface  patch  to  the  receiver,  and  the  bidirec¬ 
tional  scattering  distribution  function  (BSDF)  of  the  surface 
/, 20-22  which  is  a  dimensionless  parameter  describing 
surface- scattering  properties  that  are  invariant  in  expected 
value  to  changes  in  transmission  and  projected  area.  It  is 
noteworthy  that  /  cos  9t  cos  6r  is  simply  the  antilog  of  the 
standard  scattering  strength  of  ocean  acoustics23  S 
=  101og(/cos  0zcos  0r),  where  /  is  essentially  the  quantity 
describing  bidirectional  surface  “reflectance”  properties  that 


has  become  standard  in  radiometry  after  its  introduction  by 
the  National  Bureau  of  Standards  in  1977. 21  A  perfectly 
Lambertian  surface,  for  example,  has/=  1/77,  whereas  a  per¬ 
fectly  specular  surface  has  f=8(  —  6r)  8(  <fii  —  4>r 

+  7r)/(sin  6i  cos  0t).  The  latter  singularity  is  necessary  to  es¬ 
cape  the  implicitly  incoherent  nature  of  a  summation  over  dA 
in  Eq.  (1).  When  /  is  an  angle-independent  constant,  it  is 
related  to  the  surface  albedo24  a  by  /=  a  I  tv  and  can  be  in¬ 
terpreted  as  becoming  the  constant  /z  associated  with  Mack¬ 
enzie’s  formulation  of  seafloor  scattering  strength.23 

It  is  practical  to  recast  Eq.  (1)  in  decibels,  by  dividing 
through  by  appropriate  reference  units  and  taking  10  log  of 
both  sides.  This  homomorphic  transformation  stabilizes  the 
variance  for  optimal  pattern  recognition10,25  and  compresses 
the  dynamic  range.  The  result  is 

RdA  =  F(0i,(f>i ; dr,(pr)  +  W— TLZ  — TLr+  Ci+Cr 

+  101og(JA/Aref),  (2) 

where  the  reverberation  level  on  the  left  in  dB  re  1  /zPa  is 
balanced  by  the  scattering  distribution’s  strength  F,  the 
source  strength  W  in  dB  re  1  /zPa  @  1  m,  the  transmission 
loss  from  the  source  to  scattering  patch  TLZ ,  and  from  the 
scattering  patch  to  the  receiver  TLr ,  in  dB  re  1  m,  the  pro- 
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FIG.  13.  High-resolution  charts  of 
monostatic  and  bistatic  measured  and 
modeled  reverberation  for  S478  over 
the  region  shown  in  Fig.  2  with  upper 
and  lower  scarp  contours  overlain,  (a) 
Measured  monostatic  reverberation. 

(b)  Modeled  monostatic  reverberation. 

(c)  Measured  bistatic  reverberation. 

(d)  Modeled  bistatic  reverberation. 


jection  terms  Ct  =  101og(cos  0t),  Cr=  101og(cos  0r),  and  the 
area  term,  in  dB  re  1  m.  With  Cartesian  coordinates  for  east 
and  north  specified  by  x  and  y,  the  differential  area  becomes 
dA  =  dxdy.  In  practice,  the  elements  dx ,  dy  represent  the 
5-m  bathymetric  sample  interval  of  the  high-resolution  re¬ 
gion  of  B ' ,  over  which  the  seafloor  can  presumably  be  ap¬ 
proximated  as  planar. 

The  strength  of  the  bidirectional  scattering  distribution  F 
then  is  independent  of  the  local  seafloor  foreshortening  as 
seen  from  the  perspectives  of  the  source  and  receiver , 
whereas  the  scattering  strength  S  is  not,  so  F  is  preferred 
over  S  as  the  parameter  to  empirically  describe  the  inherent 
scattering  properties  of  the  MAR  seafloor  since  there  is  no 
unique  surface  normal  to  characterize  bathymetry  within  our 
system’s  resolution  footprint. 

To  model  reverberation,  each  of  the  spatially  varying 
terms  of  Eq.  (2),  which  excludes  W  and  101og(dA/Aref),  is 
charted  in  an  x,y  grid  over  the  high-resolution  region. 

To  compute  the  projection  terms,  range-depth  maps  of 
the  propagation  direction  of  minimum-travel-time  rays 
spreading  outward  from  a  point  source  below  a  pressure- 
release  surface  are  computed  for  a  bottomless  ocean,  given 
the  watercolumn  sound-speed  profile  shown  in  Fig.  1  of  Ref. 
2.  The  point  source  is  located  (1)  at  the  181-m  depth  of  the 
center  of  the  CORY  source  array,  for  rays  traveling  to  the 


virtual  scattering  surface,  (2)  at  the  460-m  depth  of  the  AL¬ 
LIANCE  receiving  array,  by  reciprocity,  for  rays  traveling 
from  the  virtual  scattering  surface  to  the  ALLIANCE  receiver, 
and  (3)  at  the  170-m  depth  of  the  CORY  receiving  array  for 
rays  traveling  from  the  virtual  scattering  surface  to  the  CORY 
receiver.  The  high-resolution  bathymetry  is  then  swept 
through  the  range-depth  map  to  determine  the  incident  and 
scattered  angles  at  the  bottom.  The  inner  product  of  these 
angles  with  the  local  surface  normal  yields  the  cosine  factors 
of  Eq.  (1),  and  log  transformation  yields  the  bathymetric 
projection  terms  of  Eq.  (2)  that  relate  linearly  to  reverbera¬ 
tion  level. 

Charts  of  transmission  loss  are  computed  in  a  similar 
fashion.  Range-depth  maps  of  the  monochromatic  field  pro¬ 
duced  by  a  10-element  source  array,  computed  using  the 
parabolic  equation,26  are  integrated  across  the  200-255-Hz 
frequency  band  of  the  LFM  incoherently ,  to  be  consistent 
with  actual  rough  surface  scattering.  The  resulting  incoherent 
broadband  transmission  loss,  shown  in  Fig.  15(a),  exhibits 
substantial  structure  with  a  well-defined  main  beam,  the  rem¬ 
nants  of  frequency-averaged  Lloyd  mirror  interfering  beams, 
and  sidelobes.  High-resolution  bathymetry  is  then  swept 
through  this  map  to  chart  the  incident  transmission  loss  TL* . 
A  similar  procedure  is  used  to  chart  transmission  loss  to  the 
receiver  TLr ,  except  that  a  point  source  at  the  depth  of  the 
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FIG.  14.  High-resolution  charts  of 
monostatic  and  bistatic  measured  and 
modeled  reverberation  for  S528  over 
the  region  shown  in  Fig.  2  with  upper 
and  lower  scarp  contours  overlain,  (a) 
Measured  monostatic  reverberation. 

(b)  Modeled  monostatic  reverberation. 

(c)  Measured  bistatic  reverberation. 

(d)  Modeled  bistatic  reverberation. 


given  receiving  array  is  used  to  exploit  reciprocity  from  the 
scattering  site.  This  follows  the  benign  assumption  that  mean 
transmission  loss  will  not  vary  significantly  across  the  re¬ 
ceiving  array.  The  range-depth  TLr  plotted  in  Fig.  15(b)  and 
(c)  shows  less  structure  and  leads  to  more  uniform  transmis¬ 
sion  loss  back  to  the  receiver  across  the  region  where  high- 
resolution  bathymetry  is  available,  which  lies  at  the  1/2  CZ 
vertex. 

Charts  of  the  various  projection  and  transmission  loss 
terms  of  Eq.  (2),  presented  in  Fig.  17  over  the  region  of 
high-resolution  bathymetry  for  S435,  show  substantial  spa¬ 
tial  structure  that  varies  considerably  from  monostatic  to  bi¬ 
static  geometries  when  the  bistatic  separation  is  large. 

Since  bathymetry  within  the  resolution  footprint  of  the 
MAE  bistatic  towed-array  systems  cannot  be  approximated 
as  a  planar  surface,  as  demonstrated  in  Sec.  I,  Eqs.  (1)  and 
(2)  cannot  directly  describe  measured  reverberation  and  its 
level.  The  complicated  nature  of  bathymetry  within  the  reso¬ 
lution  footprint,  however,  leads  to  a  simple  statistical  de¬ 
scription  of  the  expected  intensity  at  the  receiver  that  makes 
these  expressions  relevant.  Specifically,  the  system  resolu¬ 
tion  footprint  covers  an  annular  sector  of  roughly  50  m  by 
1000  m  at  1/2  CZ,  and  so  contains  surface  roughness  of  high 
complexity  with  respect  to  the  mean  acoustic  wavelength  of 
6.7  m.  The  total  field  received  from  the  footprint  then  be¬ 
haves  as  the  sum  of  a  large  number  of  statistically  indepen¬ 


dent  fields  scattered  from  disjoint  subregions  of  the  resolved 
seafloor  patch.27  The  expected  intensity  radiated  from  the 
resolution  patch,  proportional  to  the  variance  of  the  received 
field,  can  then  be  charted  to  the  centroid  of  the  patch  by  the 
convolution 

</(x,y|r;,rr)>=  J  J  ^G(x  —  X,  y-Y\rt,rr) 
x/(  0; .  <f>i ;  >  <£r) 

X  wltlr  cos  0i  cos  6rdX  dY ,  (3) 

where  G(x,y  | r£-  ,rr)  is  a  point-spread  function27  equal  to  the 
squared  magnitude  of  the  combined  linear  beam  pattern  and 
matched-filter  functions  that  define  the  towed-array  receiver. 
Here,  all  cosine,  transmission  loss,  and  BSDF  factors  are 
implicit  functions  of  the  spatial  coordinates  of  the  integral, 
and  the  expected  intensity  is  a  function  of  the  source  and 
receiver  locations  r* ,  rr ,  which  may  vary  in  3-D. 

The  overwhelming  contribution  to  the  convolution 
integral  typically  comes  from  the  resolution  footprint 
Af(x,y  |r?-  ,rr)  which  is  defined  in  Appendix  A.  This  obser¬ 
vation  has  been  previously  exploited  in  Refs.  28,  11,  and  2, 
and  leads  to  the  relation 
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FIG.  15.  (a)  Range-depth  slice  of 
broadband  transmission  loss  in  the 
200-255-Hz  band  modeled  incoher¬ 
ently  with  the  parabolic  equation  for 
the  Cory  10-element  source  array, 
centered  at  181  m  depth  with  2.5-m 
projector  spacing  and  MAR  vertical 
sound- speed  structure  (Ref.  2)  in  a 
bottomless  ocean,  (b)  Same  as  (a)  ex¬ 
cept  source  is  now  at  a  single  point,  at 
the  170-m  depth  of  the  Cory  receiving 
array,  (c)  Same  as  (b)  except  source  is 
a  point  at  the  460-m  depth  of  the  Al¬ 
liance  receiving  array.  The  B '  abyssal 
hill  is  typically  1/2  CZ  (—33.3  km) 
from  the  source  or  receiver.  Its  peak 
lies  at  a  depth  of  —3500  m. 


(/(x,y|r,-,rr))«  J  J  /( 6l,4>l\6r,4>r)\vlll, 

Af(x,y\ri  ,rr) 


cos  Of  cos  6rdX  dY , 


(4) 


which  is  a  very  close  approximation  to  Eq.  (3)  when  sidelobe  leakage  is  insignificant,  as  it  is  for  prominent  returns  across  a 
fixed  travel-time  ellipse.  Recasting  it  in  decibels  yields 


^  (/)(•*>  y  I r;  ,rr)s“VT+  10  log  J  J  iolF(ffi^i^r’A)-TLi-TLr+ci  +  cr]/ioc/xdY  j  _ 

\Af(x,y  |rf,rr)  / 

I - 


(5) 


The  important  issue  here  is  that  scattering  distribution’s 
strength  F  is  not  linearly  proportional  to  the  level  of  ex¬ 
pected  reverberation  unless  F  is  a  constant  over  the  resolu¬ 
tion  footprint.  Constancy  of  F  in  Ay  requires  either  it  be 
independent  of  spatial  position,  incident,  and  scattered  angle 
or  the  scattering  surface  be  a  plane  over  which  F  may  vary 
with  incident  and  scattered  angle  but  not  position.  To 
obtain  F  from  measured  reverberation  in  any  other  case  re¬ 
quires  a  deconvolution,  of  the  kind  worked  out  in  Ref.  28 
and  implemented  with  wide-area  MAE  data  in  Ref.  2.  A 
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successful  deconvolution,  however,  depends  on  a 
sufficient  set  of  independent  measurements28  as  well  as 
charting  error  that  is  insignificant  compared  to  the  dimension 
of  the  seafloor  patch  to  be  resolved.  While  both  of  these 
conditions  were  satisfied  in  the  wide-area  analysis  of  Ref.  2, 
neither  is  met  in  the  present  analysis  of  high-resolution  data. 

Given  these  considerations,  any  estimate  of  seafloor 
scattering  strength  obtained  in  the  present  analysis  will  in¬ 
evitably  be  blurred  over  the  resolution  footprint  of  the 
towed-array  system.  To  describe  this  blurring,  it  is  conve¬ 
nient  to  define 
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FIG.  16.  Spatial  distribution  of  the 
most  prominent  returns.  The  upper  5% 
of  reverberation  within  the  high- 
resolution  region  is  selected  for  each 
segment.  The  overlap  frequency  of 
these  selected  regions  is  shown  in  (a) 
for  monostatic  measured  reverbera¬ 
tion,  (b)  monostatic  modeled  rever¬ 
beration,  (c)  bistatic  measured  rever¬ 
beration,  and  (d)  bistatic  modeled 
reverberation.  Red  indicates  regions  of 
zero  intersection;  that  is,  where  rever¬ 
beration  is  only  in  the  upper  5%  for  a 
single  segment.  All  18  monostatic  and 
20  bistatic  segments  shown  in  Fig.  3 
are  used  for  this  computation. 


,rr)  =  VT+  10 log/  J  J  10(_TL'“TL''+c''+c'')/10<iZ(/r'' 

\  Afix,y\rifrr)  ) 


(6) 


as  modeled  reverberation  from  a  seafloor  that  scatters 
equally  in  all  directions.  Variations  in  modeled  reverberation 
then  arise  only  as  the  integral  of  surface  projection  and  trans¬ 
mission  loss  terms  over  Af.  This  construction  leads  to  a 
simple  linear  equation  for  the  mean  strength  of  the  scattering 
distribution  over  the  resolution  footprint 

F(x,y | r,  ,rr)  =R{I){x,y\ r, , rr) - RM(x,y | r, , rr) .  (7) 

Note  that  both  incident  and  scattered  angles  are  inte¬ 
grated  over  the  resolution  footprint.  Blurring  is  then  consid¬ 
erable,  since  these  angles  have  large  standard  deviations  over 
the  footprint,  as  shown  in  Sec.  I. 

From  Eq.  (7),  the  maximum  likelihood  estimate  for  F  is 

F{x,y\ r,  ,rr)  =  R(x,y\ri, rr)  -RM(x,y\ r ,  ,rr) ,  (8) 

namely,  the  difference  between  reverberation  measured  with 
MAE  data  R(x,y\rf  ,rr)  and  RM(x,y\ri9rr)  reverberation 
modeled  with  the  diffuse  scattering  assumption  F  =  0. 
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Under  the  assumption  that  the  receiver  measures  a  cir¬ 
cular  complex  Gaussian  random  (CCGR)  field,2,27  reverbera¬ 
tion  measured  from  a  given  surface  patch  obeys  the 
exponential-gamma  distribution.10  This  follows  from  the 
central-limit  theorem  and  the  fact  that  the  system  resolution 
footprint  is  much  larger  than  wavelength- scale  surface 
roughness.9,10,27  A  bias-corrected  estimate  for  F  is  then 

F(x,y\ri,rr)  =  R(x,y\ri,rr)-RM(x,y\ri,rr) 

+  (if/(N)  ”  In  A)  1 0  log  e ,  (9) 

following  Ref.  10  where  i//(N)  is  Euler’s  psi  function  with 
measurement  degrees  of  freedom  N  defined  in  Sec.  I. 

D.  Comparison  between  measured  and  modeled 
charts  of  high-resolution  reverberation 

Charts  of  modeled  reverberation  RM ,  obtained  from 
transmission  loss  and  surface  projection  maps  by  applying 
the  spatial  convolution  described  in  Eq.  (6),  are  presented  in 
Figs.  12- 14(b)  and  (d)  for  comparison  with  charts  of  mea¬ 
sured  reverberation  R ,  in  Figs.  12- 14(a)  and  (c). 
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FIG.  17.  Projection  and  transmission 
loss  terms  over  the  region  of  high- 
resolution  bathymetry  for  S435  from 
Eq.  (2).  (a)  Surface  projection  C, 
=  10  log(cos  0j)  of  bathymetry  from 
source  to  seafloor  using  raytrace.  (b) 
Surface  projection  Cr=  10  log(cos  0r) 
of  bathymetry  from  seafloor  to  Alli¬ 
ance  receiving  array  center  point  us¬ 
ing  raytrace.  (c)  Transmission  loss  TL ,• 
from  Cory  10-element  source  array  to 
scattering  patch  using  parabolic  equa¬ 
tion  incoherently  integrated  from  200- 
255  Hz.  (d)  Transmission  loss  TL,. 
from  seafloor  to  Alliance  receiving 
array  center  point  using  parabolic 
equation  incoherently  integrated  from 
200-255  Hz. 


Measured  and  modeled  reverberation  show  close  quali¬ 
tative  agreement  in  all  monostatic  charts.  The  most  promi¬ 
nent  returns  resolve  the  morphology  of  the  upper  and  lower 
scarps  with  reasonable  accuracy  primarily  because  the  cross¬ 
range  axis  of  the  resolution  footprint  typically  lies  along  the 
scarp  axis  and  so  does  not  bleed  into  neighboring  regions, 
such  as  the  plateau  and  trench.  Returns  from  these  neighbor¬ 
ing  regions  are  typically  about  20  to  30  decibels  lower.  Since 
this  variation  greatly  exceeds  the  expected  2.5  to  5.6  dB 
standard  deviation,  the  scarps  emerge  as  deterministic  fea¬ 
tures  of  the  data.  The  measured  and  modeled  reverberation 
charts  also  agree  on  when  the  upper  scarp  will  fall  in  the 
shadow  zone,  as  it  does  in  S528. 

The  spatial  frequency  distribution  of  the  upper  5%  of 
modeled  reverberation  for  the  monostatic  charts  is  shown  in 
Fig.  16(b).  The  modeled  distribution  bears  a  close  resem¬ 
blance  to  the  corresponding  measured  one,  with  the  weight 
of  prominent  returns  charted  to  the  steepest  portions  of  the 
upper  and  lower  scarps.  This  result  is  consistent  with  the 
high  visual  correlation  found  between  measured  and  mod¬ 
eled  monostatic  reverberation. 

The  situation  is  different  for  the  bistatic  data.  While  the 
general  character  of  measured  reverberation  shows  qualita¬ 
tive  agreement  with  the  modeled  in  all  bistatic  charts,  the 
kind  of  high  visual  correlation  found  in  the  monostatic  com¬ 
parison  only  appears  when  the  bistatic  separation  between 


source  and  receiver  azimuths  is  small  enough  for  prominent 
returns  to  faithfully  image  scarp  morphology.  More  fre¬ 
quently,  the  bistatic  separation  is  too  large  for  this,  and  el¬ 
liptic  arcs  of  both  measured  and  modeled  reverberation,  ob¬ 
liquely  centered  on  the  scarp  axes,  bleed  into  neighboring 
plateaus  over  lengths  commensurate  with  the  resolution  foot¬ 
print.  As  a  result,  the  overall  reverberation  level  in  the  vicin¬ 
ity  of  the  scarps  is  much  higher  than  in  neighboring  regions 
where  no  scarps  are  present,  as  is  often  more  clearly  seen  in 
the  wide-area  charts  of  Sec.  II  A.  This  level  fluctuates  dras¬ 
tically  along  the  scarp,  in  both  measured  and  modeled  charts, 
with  an  inner  scale  equaling  the  range  resolution  of  the  sys¬ 
tem.  The  most  prominent  modeled  returns  typically  coincide 
with  a  characteristic  nonlinearity  of  the  scarp,  such  as  a  can¬ 
yon  wall,  that  projects  area  strongly  along  the  acoustic  path 
from  source  to  scatterer  to  receiver  for  the  given  bistatic 
geometry.  The  separation  between  these  large  projections  is 
typically  much  greater  than  the  range  resolution  of  the  sys¬ 
tem,  but  shorter  than  the  cross-range  resolution.  This  is  con¬ 
sistent  with  the  greater  uniformity  in  level  found  in  mono¬ 
static  reverberation  images  of  the  scarps,  where  many 
canyons  are  blurred  together  in  cross-range. 

The  most  prominent  measured  returns  often  do  not  co¬ 
incide  with  the  most  prominent  modeled  returns  when  the 
bistatic  separation  is  large,  although  they  are  all  typically 
centered  somewhere  along  the  upper  or  lower  scarp.  To  il- 
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lustrate  this  point,  the  expected  spatial  frequency  distribution 
of  the  highest  5%  of  modeled  bistatic  reverberation  is  shown 
in  Fig.  16(d).  As  in  the  monostatic  case,  prominent  returns 
are  most  frequently  charted  to  steep  regions  of  the  scarps  in 
both  the  measured  and  modeled  distributions.  The  measured 
distribution,  however,  differs  significantly  from  the  modeled 
one  in  that  it  shows  a  relatively  high  frequency  of  prominent 
returns  charted  to  the  plateau  region  between  the  two  scarps, 
while  a  relatively  small  number  are  expected  there  according 
to  the  model.  There  are  two  plausible  explanations  for  these 
inconsistencies.  First,  the  model  predicts  variations  in  rever¬ 
beration  level  along  the  scarp  that  are  on  the  order  of  the 
5.6-dB  maximum  standard  deviation  expected  in  the  data. 
This  type  of  signal-dependent  noise,  conventionally  referred 
to  as  speckle ,  ’  ’  ’  is  then  sufficient  to  bury  the  pattern 
predicted  by  the  model.2  Second,  there  may  be  charting  er¬ 
rors  that  exceed  the  range-resolution  of  the  system.  This  is 
extremely  plausible  because  fluctuations  in  receiver  array 
heading  in  excess  of  0.1°  are  sufficient  to  cause  such  errors 
at  1/2  CZ.  Moreover,  fluctuations  of  this  order  were  often 
observed  during  reception  of  any  given  LFM  by  the  ALLI¬ 
ANCE.  The  second  explanation  is  consistent  with  an  overall 
rotation  of  predicted  reverberation  and  is  consistent  with  off- 
scarp  charting  of  prominent  returns,  but  does  not  account  for 
the  apparent  randomization  bistatic  returns  along  the  scarp  as 
does  the  presence  of  speckle  noise. 

III.  BI-AZIMUTHAL  SCATTERING  FROM  THE  B' 
ABYSSAL  HILL 

A.  Site-specific  approach 

The  development  of  Sec.  II C  makes  clear  that  it  is  typi¬ 
cally  not  possible  to  estimate  the  BSDF  of  the  B '  abyssal  hill 
directly  from  MAE  data.  Only  the  mean  BSDF,  or  equiva¬ 
lently  the  strength  of  the  mean  scattering  distribution  over  a 
resolution  footprint,  can  be  estimated.  But,  this  central  value 
is  not  necessarily  independent  of  the  method  of  measure¬ 
ment.  To  minimize  the  effects  of  the  measurement  and  to 
make  sense  of  the  bidirectional  scattering  properties  con¬ 
tained  in  reverberation  charts,  some  further  reduction  of  the 
data  is  necessary.  The  approach  behind  this  reduction  be¬ 
comes  clearer  when  the  following  four  points  are  considered. 

First,  the  seafloor  within  a  resolution  footprint  of  the 
towed-array  imaging  system  cannot  be  approximated  as  a 
planar  surface,  except  perhaps  within  the  plateau  region. 
Therefore,  we  find  it  senseless  to  plot  estimates  of  scattering 
strength  as  a  function  of  mean  incident  and  scattered  angle 
over  the  resolution  footprint,  and  instead  take  another  ap¬ 
proach. 

Second,  the  source  and  receiver  locations  for  each  trans¬ 
mission  can  be  uniquely  identified  by  the  receiver  azimuth 
which  spans  a  full  half- space  about  B ' ,  whereas  the  source 
azimuth  is  restricted  to  values  near  zero.  It  is  both  convenient 
and  natural  then  to  plot  measured  reverberation,  modeled 
reverberation  and  estimated  scattering  strength  as  a  function 
of  receiver  azimuth  with  respect  to  the  B'  axis,  with  the 
understanding  that  the  full  bi-azimuthal  description  can  be 
easily  regained  by  translation  through  Figs.  3  and  4.  For 
example,  an  estimate  of  the  mean  strength  of  the  scattering 


distribution  F(x,y  |r*  ,rr)  can  be  equivalently  written  as 
F(jc,y|ft*,ftr),  for  source  and  receiver  azimuths  I2;,  I2r 
defined  with  respect  to  the  upper  or  lower  scarp  since  there  is 
a  unique  fl, ,  12  r  pair  for  every  rt ,  rr  pair. 

Third,  the  geomorphology  of  B’  is  not  homogeneous, 
but  may  be  naturally  segmented  into  three  distinct  regions. 
These  are  the  upper  scarp,  the  lower  scarp,  and  the  plateau 
defined  in  Sec.  I.  The  results  of  Secs.  II B-D  indicate  that  the 
upper  and  lower  scarps  return  the  strongest  echoes2  because 
they  direct  the  largest  projected  areas  along  the  acoustic  path 
from  source  to  receiver.  The  plateaus,  on  the  other  hand, 
scatter  relatively  weakly  and  uniformly  since  both  incident 
and  scattered  waves  intersect  their  surfaces  at  extremely 
shallow  grazing  in  the  present  long-range  experiment.  This  is 
evident  in  both  the  measured  and  modeled  reverberation  of 
Secs.  II  B-D  when  the  cross-range  axis  of  the  system  reso¬ 
lution  footprint  is  along  the  scarp  axis.  When  the  bistatic 
separation  is  large  enough  that  the  cross-range  axis  of  the 
resolution  footprint  runs  across  the  scarp  axis,  however,  re¬ 
turns  from  the  scarps  are  often  falsely  charted  to  the  plateau 
as  described  in  Secs.  I  and  II,  and  sidelobe  leakage  from  the 
scarps  often  contaminates  true  returns  from  the  plateau 
which  are  much  weaker.  These  effects  prevent  accurate  de¬ 
termination  of  the  bi-azimuthal  scattering  properties  of  the 
plateau  when  the  bistatic  separation  is  large. 

Fourth,  measured  reverberation,  modeled  reverberation, 
and  the  estimated  strength  of  the  scattering  distribution  un¬ 
dergo  spatial  fluctuations  within  these  regions  due  to  changes 
in  geomorphology  and  the  interference  structure  of  the  re¬ 
ceived  field,  as  demonstrated  in  the  reverberation  charts  of 
the  previous  section. 

Our  approach,  in  light  of  these  points,  is  to  arrive  at  a 
single  mean  strength  for  the  bi-azimuthal  scattering  distribu¬ 
tion  across  the  upper  scarp  of  the  B '  abyssal  hill  and  a  simi¬ 
lar  mean  strength  for  the  lower  scarp.  This  is  done  by  taking 
the  average  of  F(x,y  |I2;  ,I2r)  with  respect  to  x,y  over  each 
scarp  and  plotting  the  result,  along  with  its  standard  devia¬ 
tion,  as  a  function  of  receiver  azimuth  12  r  with  the  under¬ 
standing  that  this  also  specifies  the  source  azimuth.  For  ex¬ 
ample,  the  average  of  some  function  T'Xv,};)  over  area  A  is 

<^(x,y))A  =  - J  jpr(x,y)dxdy,  (10) 

so  that  its  standard  deviation  is 

<TA{1F(x,)')}  =  sqrt((’P(x,j)2)A-('P(x,);))^).  (11) 

Since  the  average  value  (F(x,y  1 12z-  ,I2r))A  is  over  a  large 
region  A  comprising  either  the  upper  or  lower  scarp,  and  F  is 
independent  of  local  seafloor  foreshortening,  its  use  (1)  cir¬ 
cumvents  difficulties  regarding  the  use  of  scattering  strength 
to  characterize  nonplanar  and  poorly  resolved  scattering  sur¬ 
faces,  (2)  is  consistent  with  a  natural  segmentation  of  seaf¬ 
loor  geomorphology,  (3)  smooths  over  spatial  fluctuations, 
and  (4)  remains  highly  relevant  to  long-range  waterborne 
measurements  since  these  generally  entail  the  kind  of  near- 
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FIG.  18.  The  mean  reverberation  level  measured  over  (a)  the  upper  scarp 
(R(jc,y|n,-  ,flr))A  and  (b)  the  lower  scarp  (R(jc,y|n,-  ,Ilr))A  as  a  func¬ 
tion  of  receiver  azimuth  Qr  with  respective  standard  deviations 
<^{*(*,30}  and  o-Aiow{fl(x,;y)}. 

horizontal  propagation  to  and  from  the  seafloor  inherent  in 
MAE.  An  information-conserving  alternative  to  the  average 
that  has  all  of  these  properties,  but  is  far  less  wieldy,  is 
presented  in  Appendix  B. 

B.  Azimuthal  dependence  of  measured  and  modeled 
reverberation  from  the  upper  and  lower  scarps 

Curves  of  mean  reverberation  level,  measured  over 
the  upper  scarp  ,ftr)X4up  and  lower  scarp 

,flr))A  ,  are  plotted  as  a  function  of  Hr  in  Fig. 
18(a)  and  (b)  along  with  their  corresponding  standard  devia¬ 
tions  crA  {fi(x,y)},  crA  {fi(x,y)}.  These  summarize  the 

up  low 

azimuthal  dependence  of  reverberation  from  geomorphologi- 
cal  features  on  the  Bf  abyssal  hill  that  stand  out  in  remote 
acoustic  images. 

The  curve  for  the  upper  scarp  is  remarkably  flat,  while 
that  for  the  lower  scarp  shows  a  gradual  roll-off  of  at  most  5 
dB  from  the  origin  to  extreme  values  of  12  r .  Fluctuations  in 
the  mean,  however,  are  on  the  order  of  the  standard  deviation 
in  both  curves.  The  latter  ranges  from  roughly  6  dB  within 
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FIG.  19.  The  mean  reverberation  level  modeled  over  (a)  the  upper  scarp 
(RM(x,y\£li  ,Ilr))A  and  (b)  the  lower  scarp  (RM(x,y\£li  ,Or))A  as  a 
function  of  receiver  azimuth  flr  with  respective  standard  deviations 
crAjRM(x,y)}  and  o-AkJ/?M(x,.y)}. 

|I2r|<30°,  to  roughly  3  dB  outside  this  range,  so  that  the 
roll-off  in  the  lower- scarp  curve  is  barely  significant. 

The  region  of  high  standard  deviation  near  the  origin  has 
a  relatively  straightforward  explanation.  This  is  where  the 
monostatic  data  is  concentrated  and  main  beam  propagation 
from  the  source  to  the  narrow  scarps  is  most  sensitive  to 
slight  variations  in  range,  since  the  cross-range  axis  of  the 
resolution  footprint  is  along  the  scarp  axes.  This  same  sensi¬ 
tivity  partially  explains  the  complete  absence  of  a  strong 
return  from  the  upper  scarp  in  S528  monostatic  in  Fig.  14(a) 
and  (b).  Conversely,  the  standard  deviations  decrease  as  |I2r| 
increases  beyond  30°  because  the  cross-range  direction  of 
the  resolution  footprint  takes  on  such  an  oblique  angle  with 
the  scarp  axes  that  it  becomes  impossible  for  the  main  beam 
to  miss  either  scarp.  The  observed  standard  deviations  may 
also  be  partially  explained  by  minor  charting  errors  which 
would  lead  to  much  higher  fluctuations  for  small  |ftr|,  due  to 
these  same  sensitivity  issues. 

For  comparison,  curves  of  the  mean  modeled  reverbera¬ 
tion  level  over  the  upper  and  lower  scarps, 
(RM(x,y\fli,flr))A^  and  (7?M(xA|U, ,Ur))Aw  are  plotted 
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as  a  function  of  receiver  azimuth  in  Fig.  19(a)  and  (b)  along 
with  their  corresponding  standard  deviations 
frAup{/?M(x,;y)},  crA^{RM(x,y)}.  There  are  a  number  of  im¬ 
portant  similarities  and  differences  between  the  modeled  and 
measured  curves.  For  |I2r|>30°,  the  standard  deviations  are 
again  lower,  typically  near  5  dB.  Within  |I2r|<30°,  they 
increase  more  drastically  on  the  upper  scarp  than  the  lower, 
at  roughly  10  to  7  dB,  primarily  because  the  main  beam  of 
the  source  array  is  able  to  more  consistently  insonify  the 
lower  scarp  for  small  |flr|.  While  the  modeled  mean  across 
the  upper  scarp  has  a  relatively  flat  dependence  on  receiver 
azimuth,  in  accord  with  the  measured  mean,  the  modeled 
mean  across  the  lower  scarp  has  an  unambiguously  convex 
dependence,  peaking  within  |I2r|<30°  and  rolling  off  by 
roughly  10  dB  as  |ftr|  approaches  90°. 

An  explanation  for  this  behavior  requires  analysis  of  the 
projection  and  TL  terms  before  convolution.  The  mean  two- 
way  transmission  loss  (TL;(x,y  | ,ftr)  +  TLr(x,y  |ft; , 
Or))A  ,  (TLi(x,ylfli,{lr)  +  TLr(x,ylfli,flr))A,  and  mean 

up  low 

surface  projection  (C,(x,y|U; ,flr)  +  Cr(x,y\(li  ,Ur))A  , 
<c,(x,y|n,-,a  r)  +  Cr(x,y  |ftj-  over  the  upper  and 

lower  scarps  are  plotted  as  a  function  of  receiver  azimuth  for 
this  purpose  in  Figs.  20  and  21,  along  with  their  standard 
deviations. 

The  mean  surface  projection  curves  show  very  regular 
convex  behavior  with  remarkably  small  standard  deviation. 
At  the  extreme  values  of  |ftr|,  however,  levels  are  only 
about  5  dB  down  from  the  peak  value  which  lies  near  the 
origin.  Since  this  gradual  roll-off  falls  within  the  standard 
deviation  and  mean  fluctuations  of  the  two-way 
transmission-loss  curves,  which  typically  exceed  10  dB,  it  is 
lost  in  the  noise  when  the  projection  and  TL  terms  are  com¬ 
bined.  As  a  result,  mean  modeled  reverberation  exhibits  a 
relatively  constant  overall  dependence  on  receiver  azimuth, 
in  accord  with  the  measured  mean.  Extreme  fluctuations  in 
two-way  TL  arise  because  the  upper  scarp  falls  in  the 
shadow  zone  of  the  main  beam’s  refractive  path,  for  some 
segments  such  as  S528,  due  to  its  proximity  to  the  source 
array’s  conjugate  depth.  Since  the  lower  scarp  lies  much 
deeper,  both  the  standard  deviation  and  fluctuations  in  the 
mean  are  small  enough  that  a  clear  convex  dependence  arises 
in  the  mean  negative  two-way  transmission  loss  that  rein¬ 
forces  the  mean  surface  projection  trend.  This  leads  to  an 
unambiguous  central  peak  in  mean  modeled  reverberation 
level  across  the  lower  scarp  that  may  also  be  reflected  in  the 
measured  curve.  The  central  minimum  in  mean  two-way 
transmission  loss  across  the  lower  scarp  arises  because  the 
main  beam  of  the  source  array  insonifies  the  lower  scarp 
most  uniformly  when  the  propagation  path  runs  across  the 
scarp  axes.  The  relatively  small  standard  deviations  found  in 
curves  of  mean  modeled  and  measured  reverberation  level, 
compared  with  those  of  two-way  TL,  follow  from  the 
smoothing  effects  of  spatial  convolution  with  the  system 
resolution  footprint. 

The  implication  of  this  decomposition  is  that  spatial 
variations  in  transmission  loss  and  projected  area  lead  to  the 
dominant  spatial  variations  in  measured  reverberation  at 
long  range.  This  finding  for  high-resolution  bistatic  data  is 
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FIG.  20.  The  mean  two-way  transmission  loss  (a)  over  the  upper  scarp 
(TLj(x,y|flj  ,flr)  +  TLr(x,y|flj  ,Hr))Au  and  (b)  lower  scarp 
(TL/(x,y|H/  ,flr)  +  TLr(x,y|Il/  as  a  function  of  receiver  azimuth 

flr  along  with  standard  deviations. 

entirely  consistent  with  that  of  the  previous  wide-area  analy¬ 
ses  of  monostatic  data,2,11  where  the  effect  of  projected  sur¬ 
face  area  was  incorporated  implicitly  in  range-dependent 
propagation  modeling  via  the  parabolic  equation. 


C.  Bi-azimuthal  scattering  distribution  of  upper  and 
lower  scarps 

Curves  of  the  mean  strength  of  the  bi-azimuthal  scatter¬ 
ing  distribution  over  the  upper  scarp  (F(x,y\ilj  ,I2r))Aup  and 
lower  scarp  (F(x,y  |ft*  ,flr))Alow  are  plotted  as  a  function  of 
receiver  azimuth  in  Fig.  22,  along  with  their  corresponding 
standard  deviations  crA  {F(x,y)},  o"A  {F(jc,y)}.  The  mean 
strength  over  the  upper  scarp  fluctuates  about  a  constant 
value  of  roughly  —  17dB±8dB.  The  larger  standard  devia¬ 
tion  and  fluctuations  in  the  mean  for  small  |flr|  are  consis¬ 
tent  with  those  found  in  both  curves  of  measured  and  mod¬ 
eled  reverberation  level.  While  one  may  argue  that  a  very 
mild  concave  trend  is  displayed  in  the  mean  scattering 
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FIG.  21.  The  mean  surface  projection,  in  dB,  (a)  over  the  upper  scarp 
(Cj(x,y|Hj  ,£lr)  + ,nr))A  and  (b)  the  lower  scarp 
(Ci(x,y  ,ftr)  +  Cr(x,y\£li  ,nr))A  as  a  function  of  receiver  azimuth  Or 
along  with  standard  deviations. 


strength  curve  for  the  lower  scarp,  this  entire  trend  falls 
within  the  roughly  7-dB  standard  deviation  shown.  The  far 
more  important  fact  is  that  a  constant  line  may  be  drawn 
across  — 17  dB  that  falls  within  all  error  bars.  By  reciprocity, 
these  results  remain  the  same  if  the  source  and  receiver  azi¬ 
muths  are  interchanged. 

In  light  of  the  data  and  our  present  analysis,  we  conclude 
that  the  mean  bi-azimuthal  scattering  distributions  of  the  up¬ 
per  and  lower  scarps  of  the  B '  abyssal  hill  are  identical  and 
have  strengths  roughly  equal  to  the  constant  —  17dB±8dB. 
Whether  the  observed  constancy  is  due  to  innate  properties 
of  the  seafloor,  or  is  a  consequence  of  the  large  and  nonpla- 
nar  areas  integrated  by  the  resolution  footprint  of  the  towed- 
array  system,  is  not  particularly  relevant  to  the  type  of  re¬ 
mote  sensing  undertaken  in  the  present  study.  The  important 
issue  is  that  the  scarps  behave  as  if  their  innate  BSDF  is 
spatially  homogeneous  and  equal  to  the  constant  10“ 17  when 
remotely  imaged  from  waterborne  paths  at  1/2  CZ  range  or 
greater. 


FIG.  22.  The  mean  strength  of  the  bi-azimuthal  scattering  distribution  esti¬ 
mated  over  (a)  the  upper  scarp  (F(x,y\£li  ,Br))A  and  (b)  the  lower  scarp 
(F(x,y\Q,i  ,ftr))A  as  a  function  of  receiver  azimuth  along  with  the 
corresponding  standard  deviations  crA  {F(x,y)},  trA  {F(x,y)}.  The  mean 
strength  over  the  scarps  fluctuates  about  a  constant  value  of  roughly 
— 17  dB±8  dB. 

IV.  CONCLUSIONS 

High-resolution  images  of  a  typical  mid-ocean  abyssal 
hill,  called  B' ,  are  formed  from  acoustic  measurements 
made  remotely  with  a  bistatic  towed- array  system  at  1/2  con¬ 
vergence  zone  range.  For  comparison,  modeled  images  are 
generated  from  high-resolution  supporting  bathymetry, 
sampled  at  the  wavelength  scale,  with  a  combination  of  para¬ 
bolic  equation  and  ray  approximations,  under  the  assumption 
that  the  seafloor  scatters  uniformly  in  all  directions.  Narrow 
scarps,  —100  m  wide,  that  stretch  along  the  ridge  axis  for 
kilometers,  consistently  produce  the  most  prominent  features 
in  both  measured  and  modeled  images.  Measured  and  mod¬ 
eled  images  show  close  agreement  when  the  cross-range  axis 
of  the  system  resolution  footprint  is  aligned  with  the  scarp 
axis.  In  this  case,  prominent  measured  and  modeled  returns 
faithfully  image  scarp  morphology  while  neighboring  pla¬ 
teaus  produce  comparatively  negligible  returns.  This  is  true 
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even  though  the  scarps  and  plateaus  typically  share  similar 
two-way  transmission  loss  and  exactly  the  same  modeled 
scattering  strength. 

A  primary  conclusion  can  be  drawn  from  this  that  could 
not  be  drawn  in  previous  studies  that  relied  upon  lower  reso¬ 
lution  supporting  bathymetry.2,11-14  Scarps  return  the  stron¬ 
gest  echoes  because  they  project  the  largest  area  along  the 
acoustic  path  from  the  source  to  receiver.  Scattering  strength 
differences  between  the  scarps  and  plateaus,  arising  from 
their  differing  material  properties,  are  apparently  of  second¬ 
ary  importance  in  modeling  long-range  acoustic  reverbera¬ 
tion,  which  averages  over  large  inhomogeneous  seafloor  ar¬ 
eas,  since  a  good  first-order  match  between  measured  and 
modeled  reverberation  is  obtained  when  these  differences  are 
neglected. 

When  the  system’s  cross-range  axis  forms  an  oblique 
angle  with  the  scarp  axis,  prominent  returns  are  still  centered 
on  the  scarps  but  also  cross  into  neighboring  plateaus  over 
lengths  commensurate  with  the  resolution  footprint’s  cross¬ 
range  extent.  Prominent  modeled  returns  in  this  scenario 
typically  coincide  with  characteristic  anomalies  on  the 
scarps,  such  as  canyon  walls  100-200  m  wide,  that  project 
area  strongly  along  the  acoustic  path  from  source  to  receiver 
and  lie  at  spatial  minima  in  two-way  transmission  loss  to  the 
seafloor.  While  the  character  of  prominent  measured  returns 
is  similar  to  this,  a  one-to-one  correspondence  with  modeled 
returns  is  absent  at  these  oblique  angles.  The  most  likely 
explanation  for  this  discrepancy  is  that  the  expected  pattern 
along  the  scarp  is  lost  in  signal-dependent  noise  arising  from 
statistical  fluctuations  of  the  scattered  field.  Specifically,  the 
along- scarp  variations  in  reverberation  level  predicted  by  the 
model  are  on  the  order  of  the  5.6-dB  standard  deviation  ex¬ 
pected  in  the  scattered  field  by  the  central-limit  theorem. 
While  there  is  also  some  evidence  of  charting  error,  intro¬ 
duced  by  undersampled  fluctuations  in  receiving-array  orien¬ 
tation,  this  leads  to  rotational  offsets  in  the  images  as  op¬ 
posed  to  the  observed  randomization. 

This  result  can  be  cast  in  more  general  terms.  Signal- 
dependent  noise ,  commonly  referred  to  as  speckle,  is  one  of 
the  primary  factors  limiting  a  towed-array  imaging  system’s 
ability  to  remotely  resolve  seafloor  geomorphology.  Degra¬ 
dation  is  most  pronounced  when  the  expected  reverberation 
level  has  spatial  variations  that  do  not  greatly  exceed  the 
5.6-dB  standard  deviation  of  speckle  noise,  as  is  the  case 
when  the  terrain  to  be  imaged  is  relatively  uniform,  even  if 
this  means  that  it  is  uniformly  rough  as  is  the  case  along  the 
B'  scarps. 

Both  measured  and  modeled  reverberation  share  the 
same  nearly  constant  bi-azimuthal  dependence,  when  aver¬ 
aged  across  either  scarp  of  the  B'  abyssal  hill.  As  a  result, 
the  mean  bi-azimuthal  scattering  distributions  of  the  two  ma¬ 
jor  scarps  on  the  Br  abyssal  hill  are  identical  and  have 
strengths  roughly  equal  to  the  constant  —  17dB±8dB. 
Whether  or  not  this  constancy  describes  an  innate  property  of 
the  scarps  that  remains  unchanged  in  more  localized  and 
high-resolution  measurements  is  not  relevant  to  the  type  of 
remote  sensing  undertaken  in  the  present  study.  The  impor¬ 
tant  issue  here  is  that  the  scarps  behave  as  if  their  innate 
bi-azimuthal  scattering  distribution  function  is  spatially  ho¬ 


mogeneous  and  equal  to  the  constant  10“ 17  when  remotely 
imaged  from  waterborne  paths  at  1/2  +  n  CZ  range.  From 
these  ranges,  in  other  words,  the  scarps  behave  as  Lamber¬ 
tian  scatterers  with  albedo  7t/10L7. 

This  result  for  the  B '  abyssal  hill  leads  to  a  more  general 
hypothesis.  Spatial  variations  in  transmission  loss  and  seaf¬ 
loor  projected  area  cause  the  dominant  spatial  variations  in 
the  expected  value  of  long-range  ocean-basin  reverberation. 
This  hypothesis  is  appealing  because  it  means  that  long- 
range  reverberation  from  the  ocean  basin  can  be  adequately 
modeled  given  two  environmental  inputs:  (1)  the  refractive 
index  of  the  water  column,  and  (2)  bathymetry.  The  latter 
must  be  of  sufficient  resolution  to  determine  the  projected 
area  of  primary  seafloor  scatterers.  Research  is  presently  di¬ 
rected  to  test  this  hypothesis  by  analyzing  the  bi-azimuthal 
scattering  properties  of  an  inside  corner.  The  results  of  this 
research  will  be  important  because  inside  corners,  and  out¬ 
side  corners  like  the  B  ’  abyssal  hill,  comprise  the  two  major 
geological  classes  of  bathymetric  highs  in  the  world’s  mid¬ 
ocean  ridges.29 

APPENDIX  A:  THE  RESOLUTION  FOOTPRINT 

The  resolution  footprint  Af(x,y\rt  ,rr)  can  be  described 
without  loss  of  generality  in  a  simple  2-D  Cartesian  coordi¬ 
nate  system  shown  in  Fig.  Al,  where  the  receiver  rr  is  at  the 
origin  and  the  source  r,  on  the  y-axis  at  positive  a.  The 
travel-time  radius  from  the  source  is  ra  =  \jx2  +  (y  —  a)2  and 
that  from  the  receiver  is  r.  An  ellipse  of  constant  two-way 
travel  time  can  be  defined  by  holding  R  =  ra  +  r  fixed.  Points 
on  the  ellipse  are  naturally  described  by  polar  coordinates 
with  respect  to  the  receiver  location  x  =  r  cos  6  and  y 


FIG.  Al.  Elliptical  geometry  of  a  bistatic  measurement  with  a  true  feature  at 
(r,  0).  The  receiving  array  lies  along  the  x-axis  and  is  centered  at  the  origin 
so  that  the  ambiguous  azimuth  is  -  6.  If  the  receiving  array  has  azimuthal 
resolution  b(0)  =  d0,  the  system’s  cross-range  resolution  is  dS,  with  an 
ambiguity  of  arc-length  dSamh ig  as  shown. 
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=  r  sin  0  so  that  the  azimuthal  resolution  of  the  receiving  ar¬ 
ray  is  set  to  (3(0)  =  d0,  the  range  resolution  to  roughly  dr, 
and  the  cross-range  resolution  to 


/  dr 

ds  =  \l  r1  +  — - 1  dO, 


where 


dr 


\d0 


ra  cos  6 


dO  R  — asm  O' 


The  resolution  footprint  in  this  coordinate  system  is  then 


/  ra  cos  6  \ 2 

Af=dr  ds=  \!  rl+\  — - : — - 1  drdO. 


\  R.  a  sin  0 


The  cross-range  resolution  ds  then  reduces  to  rd  0  for  mono¬ 
static  geometries,  while  for  general  bistatic  measurements  in 
the  present  semicircular  experimental  geometry  about  B '  it  is 
roughly  bounded  by  ^3/2rdO^\2rdO  since  the  source- 
receiver  separation  a  never  exceeds  r. 

Minor  perturbations  in  this  elliptic  geometry  arise  when 
travel  time  range  is  converted  to  actual  spatial  range,  due  to 
the  combined  effects  of  vertical  refraction  and  bathymetric 
variation. 

APPENDIX  B:  INFORMATION-CONSERVING 
SUMMARY  OF  THE  BI-AZIMUTHAL  SCATTERING 
DISTRIBUTION  OF  THE  UPPER  AND  LOWER  SCARPS 

The  mean  reverberation  and  scattering  distribution 
strength  curves  shown  in  Figs.  18,  19,  and  22  only  summa- 
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FIG.  B2.  Histogram  of  modeled  reverberation  level  over  (a)  the  upper  scarp 
and  (b)  the  lower  scarp  as  a  function  of  receiver  azimuth  12r . 

rize  the  relevant  data  without  loss  of  information30  when  the 
samples  from  which  the  means  are  extracted  are  Gaussian 
variates.  This,  however,  is  not  the  case  with  the  present  data, 
so  an  information-conserving  summary  of  the  azimuthal 
scattering  distribution  of  the  upper  and  lower  scarps  is  pre¬ 
sented  in  Figs.  B1-B3.  Here,  histograms  of  measured  rever¬ 
beration,  modeled  reverberation,  and  scattering  strength, 
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with  samples  taken  over  the  upper  and  lower  scarps,  respec¬ 
tively,  are  given  as  a  function  of  receiver  azimuth. 
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